Power-law corrections to entanglement entropy of horizons 
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We re-examine the idea that the origin of black-hole entropy may lie in the entanglement of 
quantum fields between inside and outside of the horizon. Motivated by the observation that certain 
modes of gravitational fluctuations in a black-hole background behave as scalar fields, we compute 
the entanglement entropy of such a field, by tracing over its degrees of freedom inside a sphere. We 
show that while this entropy is proportional to the area of the sphere when the field is in its ground 
state, a correction term proportional to a fractional power of area results when the field is in a 
superposition of ground and excited states. The area law is thus recovered for large areas. Further, 
we identify location of the degrees of freedom that give rise to the above entropy. 



PACS numbers: 04.60.-m, 04.70.-s, 04.70.Dy, 03.65. Ud 

I. INTRODUCTION 

Almost four decades ago, Bekenstein realized that if 
the second law of thermodynamics is not to be violated 
in the presence of a black hole, the black-hole must pos- 
sess an entropy proportional to its horizon area [1] . The 
significance of this result became clear with Hawking's 
demonstration of black- hole thermal radiation [2] . Hawk- 
ing showed that quantum effects in the background of a 
body collapsing to a Schwarzschild black-hole of mass M, 
will lead, at late times, to a radiation of particles in all 
modes of the quantum field, with a characteristic thermal 
spectrum at a temperature equal to 
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where k is the surface gravity of the black-hole, G is the 
four-dimensional Newton's constant and is the Boltz- 
mann constant. Since the Hawking temperature fixes the 
factor of proportionality between temperature and sur- 
face gravity, one finds the Bekenstein-Hawking area law 



- -here 



(2) 



Ah is the horizon area and is the four-dimensional 
Planck length. All known black-holes in n{> 2) space- 
time dimensions satisfy the area law (AL). 

The universality of the AL has raised some important 
questions that remain unanswered: what is the dynami- 
cal mechanism that makes S'gjj a universal function, inde- 
pendent of the black-hole's past history and detailed in- 
ternal condition? Why is S-g-^ proportional to ^h? What 
is the microscopic origin of BH entropy? Are there cor- 
rections to the entropy and if so, how generic are these 
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corrections? Where are the degrees of freedom, responsi- 
ble for the entropy, located? These questions often seem 
related, which a correct theory of quantum gravity is ex- 
pected to address. 

Naturally, there has been considerable work attempt- 
ing to address some of the above questions (for recent re- 
views, see Refs. [3, 4]). Broadly, there have been two ap- 
proaches: (i) associating S^^ with fundamental degrees 
of freedom of a microscopic theory of quantum gravity 
[5] and (ii) associating S'^^ with quantum matter fields 
propagating in a fixed BH background [6-8] . 

In this work, we focus on the second approach and, in 
particular, we attribute S^^ to the entanglement of the 
quantum fields inside and outside the horizon. We show 
that it is possible to: (a) obtain generic power-law cor- 
rections to Eq. (2) which are small for large horizon radii, 
but lead to fractional powers of area for small horizons, 
(b) identify the degrees of freedom that give rise to en- 
tanglement entropy and (c) test the robustness of Eq. (2) 
and its corrections for massive quantum fields. 

Consider a quantum scalar field (in a pure state) prop- 
agating in the BH space-time. For an outside observer, 
the BH horizon provides a boundary and (s)he can only 
have information of the state restricted outside the hori- 
zon. In other words, while the full state of the field is 
pure, the state restricted outside the horizon is mixed, 
which leads to a non-zero entropy. This entropy, aka 
Von Neumann entropy, can formally be written as: 



fcgTr(plnp) 



(3) 



where p is the mixed (or reduced) density matrix. The 
above definition is the microcanonical definition of en- 
tropy will be used here. Although it is also possible to 
compute the entanglement entropy in the canonical pic- 
ture [8], its usage is restrictive due to the fact that it 
implicitly assumes positive specific heat. On the con- 
trary, however, Schwarzschild BH has a negative specific 
heat. 

About two decades ago, Bombelli et al [6] showed that 
the entanglement entropy of scalar fields is proportional 
to Ah/o?, where is the area of the boundary of the 
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region being traced over (the "horizon"), and a is the 
ultraviolet cut-off at the horizon (equivalently, the lat- 
tice spacing, when space is discretized) . Identifying a 
with the Planck length £pi, one recovers the Bekenstein- 
Hawking AL (2). These results were later reproduced by 
Srednicki [7], where by tracing over the degrees of free- 
dom inside a spherical surface of radius R, he showed 
that the entanglement entropy 



where A = AttB?. Thus, the area law can be considered 
as a consequence of the entanglement of the quantum 
fields across the horizon^. (An analytical proof of the 
area proportionality has recently been given in Ref. [10]. 
For an application of entanglement in stringy black- holes, 
see Refs. [11, 12].) 

Now, why is such a computation involving scalar fields 
in flat space-time relevant for black-holes at all? Here, we 
try to provide at least a plausible answer to this question: 
in Appendix A, we consider gravitational perturbations 
in a black-hole background, and show that for certain 
modes of these perturbations, the effective action reduces 
to that of a scalar field. Further, in Appendix B, we 
write the corresponding scalar field Hamiltonian in the 
black- hole background in Lemaitre coordinates, and show 
that for a fixed Lemaitre time, it reduces to that in a flat 
space-time. Thus for time-independent quantities such as 
entropy, scalar fields of the type considered here appear 
to play an important role. 

The computation and the area proportionality of en- 
tanglement entropy by Bombelli et al [6] and Srednicki [7] 
is based on the simplifying assumption that the quantum 
field is in the vacuum (i.e., ground) state (GS). Recently, 
two of the authors (S.D. and S.Sh.) investigated the ro- 
bustness of the entanglement entropy: area law relation 
by considering non- vacuum states [13, 14]. It was shown 
that while replacing the vacuum state by generic coher- 
ent state (GCS) or a class of squeezed states (SS) does 
not affect the AL, doing so with a class of n-particle or 
excited states (ES) results in a significant deviation from 
the AL. More specifically, if the scalar field is in a class 
of 1-particle ES, it was shown that the entropy scales as 

s ~ {^y . (5) 

where the power ji is always less than unity, and decreases 
with the increase in the number of excitations [14] (see 
the Appendix (C) for details). Thus, it was shown that 
the entanglement entropy does not always lead to AL and 



^ Although, one recovers the area law, the divergenec of the entan- 
glement entropy has been a puzzle. The reason for the divergence 
is due to the fact that the boundciry delineating the region being 
traced over is sharp [9]. 



its form crucially depends on the choice of the quantum 

state. 

Given the above results, one may draw two distinct 
conclusions: first — that entanglement entropy is not ro- 
bust — and reject it as a possible source of BH entropy. 
Second — since entanglement entropy for ES scales as a 
lower power of area — it is plausible that when a generic 
state (consisting of a superposition of GS and ES) is con- 
sidered, corrections to the Bckenstein-Hawking entropy 
will emerge. In order to determine which one is correct, 
it is imperative to investigate various generalizations of 
the scenarios considered in Refs. [6, 7, 14]. To this end, 
in this work we calculate the entanglement entropy of the 
mixed — superposition of vacuum, and 1-particle - state 
(MS). We show explicitly that the MS entanglement en- 
tropy is given by 

5 = CO (^^) [l + ci/(^)] where f{A) = (^^) , (6) 

co,ci are constants of order unity, and is a fractional 
power which depends on the amount of mixing. Thus, 
we show that, for large horizon area (^ » a^), the con- 
tribution of f{A) is negligible and the MS entanglement 
entropy asymptotically approaches the GS entropy. This 
is significantly different from the 1-particle state consid- 
ered in Ref. [14] for which the entropy always scales as 
power of area, cf. Eq.(5), the power being less than unity. 
Thus, if black- hole entropy is a consequence of quantum 
entanglement, the AL is valid for large horizons, as long 
as the quantum field is in a superposition of vacuum and 
1-particle states. 

From a physics point of view, we expect power-law cor- 
rections to Bekenstein-Hawking entropy for the following 
two reasons: (a) is a semi-classical result and is valid 
for large black- holes, i.e., when rh ^ (rt is the ra- 
dius of the horizon). It is not clear, whether the AL will 
be valid for the small BHs (r^ ^ ^pj. (b) There is no 
reason to expect that the Bekenstein-Hawking entropy 
to be the whole answer for a correct theory of quantum 
gravity. For instance, it was shown by Wald [15] that if 
one takes into account higher curvature corrections to the 
Einstein-Hilbert action, the Bekenstein-Hawking entropy 
is the leading term in a series expansion. 

As mentioned earlier, in this paper, we also identify 
the precise location of the microscopic degrees of free- 
dom (DoF) for the entanglement entropy of the superpo- 
sition of vacuum and 1-particle states [16]. We find that 
the DoF close to the horizon contribute most to the to- 
tal entropy. However, there are small contributions from 
the DoF far away from the horizon as well. These far- 
away DoF contributions are least in the case of vacuum 
state and increase as the number of excitations and/or 
the mixing weight of 1-particle state with vacuum state 
increases. Correspondingly, deviation from the AL in- 
creases as well. Thus, the corrections to the AL may, in 
a way, be attributed to the far-away DoF. 

Finally, we investigate the effect of massive fields on en- 
tanglement entropy. We show that in all cases (vacuum. 
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1-particle and superposed states), the massive scalar 
field entanglement entropy (Sm) contains an exponential 
damping factor in comparison with the massless scalar 
field entanglement entropy {So). The Gaussian fits of 
the ratio Sm/So with the mass (in appropriate units) 
show that the c;xponential factor depends explicitly on 
the mass squared and hence falls off rapidly as the mass 
is increased. However, there is not much variation of the 
fitting parameters for the different — vacuum, 1-particle 
and mixed — states, even when a fairly high amount of 
excitation is taken into account. This shows that the 
mixing proportions in the GS and ES superposition have 
little influence on the ratio Sm/So. We thus show that 
the mass overall reduces the entropy exponentially. 

The organization of this paper is as follows: In the 
next section, we briefly review the procedure of obtaining 
the entanglement entropy of massless scalar fields in flat 
space-time. In Sec. (Ill), we obtain the (reduced) density 
matrix for the scalar field which is in a superposition of 
GS and l-particlc ES. We compute the entanglement en- 
tropy numerically for such superposition and estimate the 
corrections to the BH area law. In Sec. (IV), we locate 
the scalar field degrees of freedom that arc responsible 
for the entanglement entropy for the superposition of GS 
and ES. In Sec. (V), we obtain the entanglement entropy 
for a massive scalar field. We conclude with a summary 
and open questions in Sec. (VI). In Appendix A, we dis- 
cuss the motivation for considering massless or massive 
scalar field for computing the entanglement entropy, and 
as mentioned before, we show that for certain modes of 
gravitational perturbations, the relevant action reduces 
to that of a scalar field. In Appendix B, we obtain the 
Hamiltonian of a scalar field in a general BH space-time. 
We show that this Hamiltonian in Lcmaitre coordinates, 
and at a fixed Lemaitre time, reduces to the scalar field 
Hamiltonian in flat space-time. Thus this Hamiltonian is 
relevant for the computation of time-independent quan- 
tities such as entropy. In Appendix C we briefly review 
the results obtained in the earlier works [6, 7, 14] for the 
ground state and 1-particle state. 

Before we proceed, a few comments on the notation we 
use are in order: The metric is four-dimensional with the 
signature (— , +, +, +)• We use units with kg = c = h = 1 
and set M^^ = l/(167rG'). The quantum field <^ is a 
minimally coupled scalar field. 



The Hamiltonian of massless scalar field propagating 
in flat space-time is given by Eq. (B15). In order to 
obtain the entropy, we need to discretize this Hamiltonian 
on the radial lattice with lattice spacing a. Discretizing 



the Hamiltonian such that r 



a and 



L= {N + l)a is the infrared cut-ofF^, we get 
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H = Y,Him, (7) 



Im 



where (^/^j- = ipimirj), T^im = T^im,j{rj) and 
[^im,j,'^i'm' = iSii'dmm'Sjj' ■ Up to tlic ovcrall fac- 
tor of a~^, Eq. (7) is identical to the Hamiltonian of 
A/'— coupled harmonic oscillators (HOs): 



1 ^ 1 ^ 



(8) 



where the matrix Kij represents the potential energy and 
the interaction between the oscillators = 1,...,N, 
the coordinates Xi replace the field variables ipim)- For 
the Hamiltonian (7), it is given by: 




(9) 



The last two terms denote nearest-neighbour interactions 
and originate from the derivative term in Eq. (B15). 
The most general eigen-state of the Hamiltonian (8) is a 
product of AT— HO wave functions: 



N 

V'(xi, . . . ,xn) = Y\_J^i'H„. (jc]^^ Xj^ exp 



1 



where J\fi s are the normalization constants given by 



(10) 



II. ENTANGLEMENT ENTROPY OF SCALAR 
FIELDS 



,1/4 



7ri/4 ' 



(11) 



In this section, we briefly review the procedure of ob- 
taining entanglement entropy for scalar fields propagat- 
ing in flat space-time. The motivation for considering 
scalar fields for the entanglement entropy computations 
is discussed in Appendix A. The relevance of the scalar 
field Hamiltonian in flat space-time for computing en- 
tropy in a black-hole space-time is discussed in Appendix 
B. 



X = Ux, {U'^U = In), x^ = {xi,...,xn), x"^ = 
(Si; • ■ • jXat)) Kd = UKU^ is a diagonal matrix with 
elements koi, and Ui{i = 1. . .N) are the indices of the 



^ In discretizing the terms containing the derivatives, one usually 
adopts the middle-point prescription, i.e., the derivative of the 
form f{x)dx{g{x)] is replaced by fj+i/2[93+i - QjMa- 
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Hcrmitc polynomials {Hu)- The frequencies are ordered 
such that koi> koj for i > j. 

Defining the N x N matrix Q = U'^K^^U, such that 
|0| = lii'ol-^/^, and tracing over first n of the A'' oscilla- 
tors, one obtains the reduced density matrix: 

/n 
Yldxi -ipixi, ...,Xn;t) ■4'*ixi,. . .,Xn;t') 

JJdxiexp -^-^ H-^*^'^- 
i=i L -I i=i 

AT 

U^j-H.^i'^m^^ (12) 



X exp 



x'^flx' 



3=1 



where we now denote: x = (xi, . . . , Xn;ti, . . . , tN-i) = 



(xi 



,Xn;t), with t = ti, 



tN-r 



J 



1, . . . ,{N — n). It is easy to check that p, implying 
that p is mixed, i.e, although the full state is pure, the 
state obtained by integrating over n-HO is mixed. Sub- 
stituting the reduced density matrix (12) into the formal 
expression (3) will yield a non-zero (positive) entangle- 
ment entropy. 

It is not possible to obtain a closed form expression for 
the density matrix for an arbitrary state (10). However, 
in the cases where all the HOs are in their GS [7], or in 
the GCS or in a class of SS [14], all of which are mini- 
mum uncertainty states, closed form analytic expressions 
of p{t; t'), and hence of the entropy, can be evaluated ex- 
actly and shown to follow the BH AL. For the first ES, 
not a minimum uncertainty state, the entropy computed 
numerically [14] is found not to obey the AL. 

In the following section, we obtain the entanglement 
entropy for the superposition of GS and ES. [For the sake 
of completeness, we have briefly discussed the entangle- 
ment entropy for the ground and first-excited states in 
the Appendix C.] 



III. ENTANGLEMENT ENTROPY FOR A 
SUPERPOSITION OF GS AND ES 

In this section, we obtain the entanglement entropy for 
the superposition of ground and excited states. (In the 
following, we denote all relevant quantities such as the 
wave function, density matrix etc by the symbol/suffix 
for GS and by 1 for the first ES.) 

The discretized scalar field wave function ^/i in a MS is 
a linear superposition of the N-HO GS wave function ipo, 
Eq. (CI) [Appendix C], and N-HO (1-particle) ES wave 
function tpi [corresponding to one HO in the ES, while 
the rest N-lin their GS, Eq. (C12)], i. e., 

ip{x;t) = [cotpo{x;t) + a tlji{x;t)] (13) 

where x = {xi, ■ ■ ■ ,Xn} ; and as before tj = Xn+j {j = 
1, - ■ ■ ,N - n) ; t = {ti, ■■■ , ijv-n} = {xn+i, ■■■ , xn}- 



Normalization of ip requires CQ + cf = 1. Here we assume 
that Co and Ci are real constants. 

Referring to the Appendix C and using Eq. (C12), we 
can write, 



where, 



i){x;t) = [co + ci f{x;t)]Mx;t), (14) 



f{x;t) = V2a'^K]^^Ux = y^x, (15) 



a being the expansion coefiicient defined in Eq. (C14). 
The N dimensional column vector y is given by 



y = V2U^K]/\ 



VA 

Vb 



(16) 



yA and ys are n- and (A^ — n)-dimensional column vec- 
tors, respectively. 

The density matrix is a sum of three terms: 



« n 

Pit;t') = J l[dxi ^{x;t)^*{x;t') 



(17) 



= cl po{t;t') + c\ pi{t;t') + coci P2{t;t') 

where po{t\t') is the GS density matrix (C4), pi{t-,t') is 
the ES density matrix (C15). It is easy to see that, one 
can make the following identifications of the matrix A 
and it's components; and the constant k (see Appendix 
C), with the column vector y and it's components as 



A - T / Ayi As 

^ = --yy = ' Ac 



= ; As = iy^ys ; Ac = ^ysys , 

K = Tt{AaA-^) = ^ylA-^yA . 



(18) 



P2 is the cross-term in the total density matrix p, Eq. 
(17), due to the mixing of GS and ES and can be evalu- 
ated as follows: 



/n 
[] dXi [f{x; t) + fix; t')] Mx; t) Vo(*; t' 



[VB-pf {t + t')po{t;t'), 



where 



p = B^A-^yA , 



(19) 

(20) 



is an [N — n)-dimensional column vector"^. 

Using Eqs. (C4), (C15) and (19), the complete MS 
density matrix can be written as 

p{t; t') = [cl + cIk {1 + u{t; t')} + coCiv{t] t')] po{t\ t') 

(21) 



For definitions of matrices A, B etc. see the Appendix C [Eq. 
(C2)]. 
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where the functions u and v are defined by 



and 



u{t;t') = - t^^^t+^t'^^^t' ^ ^T^^^, 
v{t;t') = {vB-pf (t + f). 
Let us define 



(22) 



F{t;t') = 1 + Kiw{t;t') + K2v{t;t') + -j^^{t;t') (23) 



where 



w{t;t') = - t^^yt+f^^yt' ^ ^T^^^^, 



Ac 



( Ac 
A^/ = A^ + 2k;o A/3 - — 



(24) 



and 



Ko = — ; «;! = — ; K2 = — , « = cq + cjk . (25) 

K K 

A^' and A^' arc {N — n) x (N—n) matrices, and constants 
{ko,ki,K2) describe the amount of mixing between the 
GS and ES. 

With these definitions, the density matrix (21) can be 
rewritten as 



p{t;t') = K F{t;t') po{t;t') 



(26) 



As for the ES, hero too the prc-factor F{t; t') of the Gaus- 
sian po{t;t') cannot be factorized into {N — n) two HO 
density matrices. However, as discussed in Appendix C 2, 
if the vector t'" is outside the maximum t^ax^ given by 
Eq. (C19) corresponding to the Str Hmits, the argument 
of po{t] t') is neghgible. Therefore, if the conditions (C18) 
as well as the conditions 

are satisfied, then we can approximate the pre-factor 
F{t; t') as 

F{t\t') K. exp[Ki «;(t;i') + K2v{t\t')], (28) 

where we have kept terms up to quadratic order in 
[Note that, v{t; t') is only linear in t, t' whereas w{t; t') is 
quadratic in t,t' ] 

Using Eq. (C4) for po{t;t') we can now write the (ap- 
proximated) MS density matrix as 



pit;t') = ki 



\n\ 



7r^-»|A| 



exp[z{t;t') +k2v{t;t')] (29) 



where 



z{t; t') = - ^V^+^'^y^^ ^ ^T^,^, (30) 



= 13 + kxKjS' =13 + kiA0 - 2koki A^ - 



7' = 7 + KiA^ 



7 + KiA^ + 2koki [A/j — ) (31) 



Ac 

K 

Ac 



are {N — n) x {N — n) matrices. /3' is symmetric while 
7' is not necessarily symmetric. 

Let us make the following transformations on the set 
of {N — n) variables t = {a;„+i,--- ,xn} and t' = 



t 



t + s ; t' ^ t' + s 



(32) 



where s = {si, ■ ■ ■ ,SN-n} is a set of {N — n) constant 
values. The density matrix (29) reduces to 



p{t;t') = Af exp 



where the normalization constant Af is given by 
Af = k\ 



\n\ 



rN-n 



exp 



(33) 



(34) 



The (A'' — n) -dimensional constant column vector s is 
determined from the equation 



(35) 



It is easy to check that for either cq = or ci = 0, 
the constant ^2 vanishes, whence s = 0, and the density 
matrix (33) reduces either to that of pure GS [7] (for 
Co = 1, ci = 0, whence /?' = /3, 7' = 7) or that of ES [14] 
(for Co = 0, ci = 1, whence (3' = (3 + Afj, 7' = 7 + A,^). In 
general, when both cq and ci are non-vanishing, then un- 
der the shifts /3 — > /3', 7 — > 7' (where and 7' are given 
by Eqs. (31)) the MS density matrix (33) is of the same 
form as the GS density matrix (C4), up to a normaliza- 
tion factor given above. Such a normalization constant 
does not affect the entropy computation. Therefore we 
can use the same steps as for GS [Eqs. (C5) - (C9), with 
the replacements /? — > /3',7 — > 7'] to calculate the total 
MS entropy. 

The rest of the analysis in this section is similar to 

that of Rcf. [14]. We compute the entanglement entropy 
numerically (using MATLAB) in each of the cases: 

(i) GS (co = l,ci=0), 

(ii) ES (co = 0,ci = l), 

(iii) an equal mixing (MS^g) of ES with GS (cq = Ci = 
I/V2), and 

(iv) a high mixing {MShz) of ES with GS (cq = 1/2, ci = 
73/2). 
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In S vs In (R/a) plots for GS, ES & MS(Eq/Hi) 




4.8 5 5.2 4.8 5 5.2 4.8 5 5.2 

In (R/a) In (R/a) In (R/a) 

FIG. 1: Plots of logarithm of GS, ES and MS (Eq/Hi) entropies versus ln(_R/a), where R — a{n+ 1/2) radius of the hypothetical 
sphere (horizon), for A'' = 300, n = 100 — 200 and o — 30, 40, 50 (in the ES and MS cases). The numerical precision is 0.01%. 



The computations have been done with a precision^ of 
0.01% for the set of values: N = 300, n = 100 - 200 
and o = 30, 40, 50, o being the last non-zero columns of 
the vector . The conditions (C18) as well as (27) are 
satisfied for these values of the parameters. 

The expectation value of energy, £, for MS can be ex- 
pressed as 



2 ^ 

o ^ — ' 

i=N-o+l 



.1/2 



(36) 



where £o = | X]i=i the (zero-point) GS energy. 

Now the excess of energy over the zero point energy is 
given by the second term in the above equation (36). As 
the value of ci is between and 1 and since the /cd^'s are 
in ascending order (kui > kj^j for i > j), the fractional 
change in energy (£ — £o)/£ is at most about ~ 5% (cor- 
responding to the extreme situation ci = 1, i.e., ES), for 
iV = 300 and o ~ 50. Moreover, since there are o num- 
ber of terms in the sum in the second term of Eq.(36), 
the excitation energy {£ — Sq) is of the order unity (in 
units of 1/a, where a is the lattice spacing). Therefore if 
a is chosen to be of the order of Planck length, then the 
above energy is of the order of Planck energy. The mass 
of a semi-classical BH, on the other hand, is much larger 
than the Planck mass. Hence, one may safely neglect 
the back-reaction of the scalar field on the background 
geometry. 

In Fig. (1), we have plotted the logarithm of the en- 
tropy S versus ln(i?/a) = ln(n -I- 1/2), for different val- 
ues of the excitation (o = 30,40,50), for GS, ES and 
MS (Eq/Hi). For GS, the plot is very nearly the same 
as the numerical straight line fit obtained in Ref. [7], 



* The computations here are one order of magnitude more accurate 
than those in Ref. [14] 



S = 0.3(i?/a)2 with iV = 60 lattice points. For the MS 
(Eq/Hi) cases, as well as for ES, the plots are nearly 
linear for different values of the excitations o = 30, 40, 50 
and appear to coincide with the plot for GS for large areas 
{A = AnE? ^ a^). Numerical straight line fittings of the 
logarithm of the ES entropy, Ses, with ln(i?/a) shown in 
[14] revealed that for smaller areas Ses ^ where /i is 
always < 1 and decreases as the number of excitation o 
increases. To have a closer look at the behaviour of MS 
entropy, Sms, (foi' both equal and high mixings) and the 
ES entropy Ses with respect to the GS entropy, Sgs, we 
have plotted in Fig. (2) the ratios Sms / Sgs, Ses/ Sgs 
and the inverse ratios Sgs/SmSi Sgs/SeSi versus the 
area A. For the range of excitations (o — 30, 50), all the 
ratios tend to unity as the area increases. Thus the gen- 
eral criterion of 'asymptotic equivalence' [17] is fulfilled. 



lim , 



SxsiA) 



= 1 ; 1™^^= 



Sgs (A) 



= 1 



(37) 



Sgs{A) ^ ' '""^-^ SxsiA) 

where XS = MS (Eq or Hi) or ES. In other words, the 
MS (Eq/Hi) and the ES entropies coincide asymptoti- 
cally with the GS entropy. However, as is evident from 
Fig. (2), the MS(Eq) entropy is closer to the GS en- 
tropy for large A, than the MS (Hi) entropy and the ES 
entropy, the latter being the farthest. This implies that 
the asymptotic behaviour is strongly influenced by the 
relative weight ci of the mixing of ES with GS — the 
smaller the value of ci the sharper is the asymptote. 

In order to make things more transparent we have plot- 
ted in Fig. (3) the best fit ratios of the MS entropies (for 
equal and high mixings, with o = 30, 40, 50) to the GS 
entropy using a simple formula: 



Sms 
Sgs 



CTO 



(38) 



The fitting parameters ctq , cr and v are shown in Table 
1. The parameter ctq is very close to unity for all values 
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1.15 



Asymptotic behaviour for o = 30 
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Asymptotic beliaviour for o = 50 
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FIG. 2: Plots of ratios of GS and MS (Eq/Hi) or ES entropies and their reciprocals versus the area A (in units of a^, a being 
the lattice spacing) for o = 30, 50. The plots show the asymptotic nature of the MS and ES entropies with respect to the GS 
entropy. The curves on the upper half (above 1) show the variation of Sxs/Sas with A, where XS stands for MS(Eq/Hi) or 
ES, while the lower curves show the variation of Sgs/Sxs with A. 



Relative MS(Eq) entropy 



Relative MS(Hi) entropy 




X 10 



X 10 



FIG. 3: Best fit plots (solid lines) of the relative mixed state entropies {Sms / Sgs) for equal and high mixings versus the area 
A (in units of a^), for o = 30, 40, 50. The corresponding data are shown by asterisks. 



Fitting 
Parameters 


For MSb, 


For MSh,; 


o = 30 


o = 40 


o = 50 


= 30 


o = 40 


= 50 


CTO 


1.001 


1.002 


1.003 


1.001 


1.004 


1.006 


a 


1738 


4288 


8039 


2956 


7652 


14120 


V 


1.180 


1.210 


1.225 


1.141 


1.178 


1.192 



TABLE L Values of the parameters of the fit Sms I Sgs =0-0+0- {Aja^) " for both MS(Eq) and MS(Hi) cases with the 
amounts of excitation o = 30, 40, 50. 



of o = 30,40,50, for both MS(Eq) and MS(Hi) cases. 
Ho-wever, there is a slight increase in cto as o increases or 
for greater relative weight c\ of mixing of ES with GS 
for a particular o (i.e., co is greater in the MS (Hi) case 
than in the MS(Eq) case for fixed o). Neglecting this 
variation in (Tq and noting that the GS entropy can be 
written as Sqs = no(^/a'^), where no is a constant, we 



can approximately express: 

Sms — Sgs + ( ^ ) ' ("^^^ 

where a — riQa. As the value of exponent (1 — i^) lies 
between and —1 for both equal and high mixings (see 
Table 1) the second term in the above Eq. (39) may 
be regarded as a power law correction to the AL, result- 
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FIG. 4: Plot of the distribution of entropy per partial wave [(21 + l)^;] in the cases of GS, ES and MS (Eq/Hi), for A'^ = 
300, n = 100 and o = 30, 40, 50. 



ing from entanglement, when the wave- function of the 
field is chosen to be a superposition of GS and ES. It 
is important to note that the correction term falls off 
rapidly with A (due to the negative exponent) and in the 
semi-classical limit {A ^ 1) the AL is hence recovered. 
This lends further credence to entanglement as a possible 
source of black-hole entropy. The correction term is more 
significant for higher excitations o or greater ES-GS mix- 
ing proportion ci. This is evident from Table 1, which 
shows that the parameter a (and hence a) increases and 
the parameter i' (and hence the negative exponent |1 — 
decreases with the increase in o (fixed ci) or the increase 
in ci (fixed o). 

Fig. (4) shows the variation of {21 + l)Si with I, in the 
cases of GS, ES and MS (Eq/Hi) for a fixed n{= 100) 
and a set of increasing values of o. For the GS, there 
is a peak at Z = (s-wave), followed by another one at 
I « 40 due to the degeneracy factor {21 + 1). The first 
peak shifts to a value ? > for the ES, and the shift is 
greater as o is increased. There is, however, no second 
peak in this case, although there seems to be an increase 
towards higher values of I. Thus, higher partial waves are 
seen to get excited with greater excitations. In each of 
the MS cases, there is a trace of the first peak at / = as 
for GS, however the amplitude of that peak is very small 
compared to the second peak which appears between I ~ 
10 — 30 depending on the value of the excitations o. As 
in the case of ES, the second peaks for MS (Eq/Hi) are 
higher and far away from / = for increasing values of o. 
However, relative to the ES case, there is a broadening of 
the half- width of the peaks for MS, though not as broad 
as that of the second peak for GS. Thus, as expected. 



the {21 + l)Si vs I curves for MS show features that are 
intermediate between those for GS and ES. 



IV. LOCATION OF THE DEGREES OF 
FREEDOM 

Let us now examine closely the expression for the in- 
teraction matrix Kij, Eq. (9), for the system of N HQs. 
The last two terms which signify the nearest-neighbour 
(NN) interaction between the oscillators, are solely re- 
sponsible for the entanglement entropy of black-holes, 
i.e., Sbh = if these two terms are set to zero. Let 
us, however, consider the situation where the NN inter- 
actions, and hence the off-diagonal elements of Kij, are 
set to zero (by hand) everywhere except in a 'window', 
such that the indices i,j runs from q — s to q + s, where 
s < q. Thus the interaction region is restricted to a 
width of d = 2s + 1 radial lattice points. Now, choosing 
the position of the center of the window q to vary be- 
tween and a value qmax > n, we allow the window to 
move rigidly across from the origin to a point outside the 
horizon. Fig. (5) shows the variation of the percentage 
contribution of the entropy for a fixed window size of 5 
lattice points (c? = 5, s = 2), i.e., 

pc{q) = X 100 (40) 

Jtot 

as a function of q for fixed values N = 300, n — 100 in 
each of the cases GS and ES, MS (Eq/Hi) with o = 30, 50. 
Here Stot is the total entropy with all the NN interactions 
present, i.e., i,i running from to N. 
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GS for n = 100 



ES for n = 100 




FIG. 5: Plots of the percentage contribution pc{q) to the total entropy as a function of window position q, for a window size 
d = 5 and fixed iV = 300, n = 100, in each of cases of GS, ES and MS (Eq/Hi). For ES and MS (Eq/Hi) the sohd curve is for 
o = 30 whereas the broken curve is for o = 50. 



In all the cases of GS, ES and MS (Eq/Hi) the first 
observation is that pc{q) — when q is far away from n. 
There is no contribution to the total entanglement en- 
tropy if interaction window does not include the horizon. 
For values of q very close to n there are significant contri- 
butions to Stot and in the case of GS, pc{q) peaks exactly 
at g = n. For ES and MS, however, the peaks tends to 
shift towards a value q > n, its exact position depending 
on the amount of excitation o. Moreover, the amplitudes 
of the peaks gradually diminish as the value of o and/or 
the mixing weight ci increases. Thus, we infer that: 

• The contribution to the total entropy is more from 
the DoF that are in the vicinity (inside or outside) 
of the horizon, rather than far from it. 

• The contributions, however small, from the DoF 
far away from the horizon are more for MS (Eq/Hi) 
and ES, compared to the GS. In other words, the 
contributions from the far away DoF increases with 
increase in the number of excitations and amount 
of mixing of ES with GS. 

Further investigations have been carried out recently in 
Ref. [16] to check the effects of the far-away DoF on the 
total entropy, by keeping fixed the center of the window 
at the horizon, i.e., q = n, while increasing the window 
width d from to n. It is found that for GS about 85% 
of the total entropy is obtained within a width of just 
one lattice spacing, and within a width of d = 3 almost 
the entire GS entropy is recovered. Thus most of the 
GS entropy comes from the DoF very close to the hori- 
zon and a small part (about 15%), has its origin deeper 



inside. For ES, however, the corresponding figures are 
about 60% {d = 5), and the total ES entropy is recov- 
ered when d is as much as 15 — 20, depending on the 
number of excitations o = 30 — 50. Thus the far-away 
DoF contribute more to the entropy for the ES. This, in 
turn, may be looked upon as follows: larger the deviation 
from the area law, larger is the contribution to the total 
entropy from the DoF that are far away from the horizon. 
The situation is intermediate for the MS (which itself in- 
terpolates between the GS and ES) : This is evident from 
Fig. (5) [and also from Fig. (1)] where unlike the curves 
for ES, those for the MS cases do not show much devi- 
ations from the curve for GS, even for high excitations 
o. 



V. ENTANGLEMENT ENTROPY OF MASSIVE 
SCALAR FIELD 

As shown in Appendix A, the equation of motion for 
metric perturbations in a general space-time with a cos- 
mological constant |A| coincides with that of a test mas- 
sive scalar field propagating in the background metric. In 
all our earlier analysis, we had set, for simplicity, |A| — 0. 
In this section, we obtain the entanglement entropy for 
the massive scalar field. 

The action for the massive scalar field (mass m) prop- 
agating in the background space-time g^^ is 

S^-\ j d^x^ [g^''' d^v d,^ + m2(p2] (4^) 
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S /S„ for GS (n = 1 00) S /S„ for ES (n = 1 00, o = 50) 

m ^ ' m ^ ' 




FIG. 6: Best fit plots of the relative variation of the total entropy Sm for a massive scalar field (in units of the total entropy 
So corresponding to a massless scalar field) with the mass m times the lattice spacing a, for fixed n — 100, o — 30, in each of 
cases of GS, ES and MS (Eq/Hi). The corresponding data are shown by asterisks. The fits show an exponential damping of 
the ratio Sm/So with mass. 



Repeating the steps described in Appendix B for the 
massive scalar will lead to massive, free field Hamilto- 
nian (B15). Discretizing the resulting Hamiltonian, as 
described in Sec. (II) will lead to the A^— coupled HO 
Hamiltonian, with the interaction matrix -R'l™'' given by: 



K. 



(m) 



Kij + (ma) 



(42) 



where Kij is the interaction matrix, given by Eq.(9), for 
massless scalar field. Following the steps discussed in 
sections (II) and (III), we can compute the entanglement 
entropy for the massive field (Sm)- 

In Fig. (6), we have plotted Sm/So [where, Sq is 
the entropy for the massless scalar] for the GS, ES and 
MS(Eq/Hi) for N = 300, n = 100, o = 50. The Gaussian 
fits indicate an exponential fall-off of Sm with respect to 
Sq as the mass increases: 



Sm — Sq exp — Qfi (ma + 02) 



(43) 



where ai, (Xi and A are the fitting parameters. Depending 
on the state (GS, ES or MS), the parameter a\ varies 
between 2.77 and 3.24, is between 0.077 and 0.001 and 
the power A is close to 2. Thus approximately Sm/So 
scales as e~™ ° . There is, however, a small variation 
in the power A for the different cases. The exponential 
damping is strongest for GS, and gradually slows down 
as more and more ES oscillators are mixed with GS, the 
damping is slowest for the ES case. 

Although Sm scales as So times a mass-dependent 
exponential term, the fitting parameters a\,ai and A 



change very little for the different cases GS, MS (Eq/Hi) 
and ES, even for a fairly high amount of excitation 
o — 50. As such, for a fixed mass m, the variations In S 
versus In (i?/a) for all the cases remain almost the same 
as those for massless scalar field, cf. Fig. (1). The anal- 
ysis and inferences of the previous sections for massless 
scalar go through for the massive scalar field, resulting 
in correction terms obtained before. 



VI. CONCLUSIONS 

In this work, we have obtained power-law corrections 
to entanglement entropy, which may be relevant for the 
entropy of BH. Indeed, as shown in Appendices A and 
B, certain modes of gravitational perturbations in black- 
hole space-times behave as minimally coupled scalar 
fields. Also for computations of time-independent quan- 
tities done at a fixed value of Lemaitre time, it suffices to 
consider an effective flat space Hamiltonian. Extending 
the analysis of the earlier work [14] , we have shown that 
for small black-hole areas the area law is violated not 
only when the oscillator modes that represent the scalar 
DoF are in ES, but also when they are in a linear super- 
position of GS and ES. We found that the corrections to 
the AL become increasingly significant as the proportion 
of ES in the superposed state increases. Conversely, for 
large horizon areas, these corrections are relatively small 
and the AL is recovered. 

It is interesting to compare and contrast the power-law 
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corrections obtained here to those derived in the case of 
higher-derivative gravity [15]. The power-law corrections 
to the Bekenstein-Hawking entropy derived in both — 
entanglement and higher derivative gravity - these ap- 
proaches have same features. For instance, it was shown 
that the entropy of 5-dimensional Boulware-Deser black- 
hole [18] is given by 



will it be consistent with the first law of BH thermody- 
namics? (iii) Are the second and third laws of thermody- 
namics valid for this entropy? (iv) Can the entanglement 
of scalar fields help us to understand the evolution or dy- 
namics of BHs and the information loss problem? We 
hope to report on these in future. 



S = 



1 + 



A2/3 



c = constant . 



(44) 
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As in Eq. (39) the above entropy is proportional to area 
for large horizon radius, however it strongly deviates in 
the small horizon limit. It is important to note that 
the corrections to the black-hole entropy are generic and 
valid even for black-holes in General relativity without 
any higher curvature terms"''. It is interesting to investi- 
gate the relation between the entanglement entropy with 
that of Noether charge approach [19] . 

The location of the DoF that give rise to the entan- 
glement entropy has important implications as far as the 
corrections to the AL is concerned. It is found that for 
GS, ES and MS, the contributions to the total entropy 
from the DoF that are nearest to the horizon are maxi- 
mum. However, there are small contributions from DoF 
that are far away from the horizon, which also need to 
be taken into account in order for the AL to emerge. 
These contributions arc least in the case of GS and grad- 
ually increase as the proportion of mixing of ES with GS 
and/or the amount of excitation increases. Correspond- 
ingly, there are increasing deviations from the AL. Thus 
one is led to conclude that the AL is intimately linked 
with near horizon DoF. 

We have also shown that the mass of the scalar field 
does not have much influence on the corrections to the 
AL. The total entropy for the massive field scak;s as that 
of the massless field times a mass-dependent exponen- 
tially damping term that varies very slowly with the mix- 
ing proportion and the amount of excitation which are 
key to produce the AL corrections. 

Open problems in the context of entanglement entropy 
include: (i) The proportionality constant in the relation 
S = 0.3{R/af for GS obtained in ref.[7] differs from the 
1/4 in the Bekenstein-Hawking relation [Eq.(2)]. This 
discrepancy persists for MS and ES. A probable reason 
behind this mismatch is the dependence of the prc-factor 
on the type of the discretization scheme. For example, 
another discretization scheme, resulting in the NN in- 
teractions between four or more immediate neighbours, 
woiild result in a different prc-factor. Is it then at all 
possible to obtain the Bekenstein-Hawking value? (ii) 
Can a temperature emerge in the entanglement entropy 
scenario, and if so, then along with the current entropy, 
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Note added: Recently, in Rcf. [20], the authors have 
calculated the sub-leading power-law corrections to the 
Bekenstein-Hawking entropy using the canonical ensem- 
ble (aka Brick- wall) approach [8]. The results reported 
there agree with the numerical results derived in this 
work. 



APPENDIX A: 



WHY CONSIDER SCALAR 
FIELDS? 



In this appendix, we briefly discuss the motivation 
for obtaining the entanglement entropy of a scalar field. 
First, we obtain the equation of motion of the metric per- 
turbations for a general space-time and then, as a special 
case, discuss the equation of motion of the perturbations 
in asymptotically flat spherically symmetric space-times. 

Consider the Einstein-Hilbert action with a positive 
cosmological constant ([A]): 



(Al) 



Let 

— 9iiv + ^iiv expand the action keeping only 
the parts quadratic in h^v, we get [21] 



|A1 



where 



-(V^/i« + V./i«-V«V) 



(A2) 

(A3) 
(A4) 



The above action is invariant under the infinitesimal 
gauge transformation /i^j^ — > /ipi/-FV(^^^) when the back- 
ground metric satisfies the vacuum Einstein's equa- 
tion with |A|. We can remove the gauge arbitrariness by 
imposing the harmonic gauge condition dji'^" = [21]. 
Assuming /i^^ to be small, we can keep only the first 

as the numerical error we obtain is much larger than ln(n-|- 1/2). derivatives of hn,u- With these twO Conditions, the action 



^ In this context, it should be mentioned that it is not possible to 
check for logarithmic corrections to the entropy in our analysis. 
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(A2) reduces to (for more details see Ref. [22], pgs. 330- 
332) 



M2 

_l 

2 



?! [v«vv«/i''" + |A|v/i''"]- 

(A5) 

The above action corresponds to the massive spin-2 h^^i, 
field propogating in the background metric g^^j/ where the 
cosmological constant appears as mass-term. In the weak 
field limit — when the gravitational field is weak like in 
the case of regions close to the black-hole horizons — /i'^'' 
can be approximated as a plane-wave perturbation with 
a particular frequency, i. e. h^^, = M~^e^v'p{x^) [where 
e^v is the constant polarization tensor] , the above action 
can be written as 

S.A9,h) = Jd'x^l [9„^a>+ |A|^2] _ (A6) 

which is the action for the massive scalar field propagat- 
ing in the background metric g^v 

Moreover, for four-dimensional spherically symmetric 
space-times, the metric perturbations are of two kinds 
— axial and polar [23-25]. The equations of motion of 
both these perturbations are scalar in nature and are re- 
lated to each other by a unitary transformation [23] . The 
equation of motion of the axial perturbations is nothing 
but the equation of motion of a test, massless scalar field 
propagating in the black-hole background: 



1 



'-9 



5^ = 0. (A7) 



Hence, by computing the entanglement entropy of the 
scalar fields, we obtain the entropy of a class of metric 
perturbations of the background space-time. Of course, 
a generic perturbation being a superposition of plane 
wave modes, and entanglement entropy being a non- 
linear function of the wave-function, we do not claim that 
such a computation would account for the entropy of all 
perturbations. Nevertheless, it is expected to shed im- 
portant light on the role of entanglement in the AL. 

In most part of this work, we calculate the entropy of 
the massless scalar field (i.e., setting |A| = 0). In Sec. 
(V) , we obtain the entanglement entropy of the massive 
field, corresponding to |A| ^ 0. 



APPENDIX B: HAMILTONIAN OF SCALAR 
FIELDS IN BLACK-HOLE SPACE-TIMES 

In this appendix, we obtain the Hamiltonian of the 
massless scalar field propagating in a general static spher- 
ically symmetric space-time. We show that for a partic- 
ular choice of time slicing such a Hamiltonian reduces to 
the Hamiltonian of a scalar field in flat space-time. 

In Ref. [16], two of the authors (SD and SSh) showed 
that in a fixed Lemaitre time coordinate, the Hamiltonian 



of the scalar field propagating in Schwarzschild space- 
time reduces to the scalar field Hamiltonian in flat space- 
time. In this appendix, we extend the analysis for any 
non-degenerate static spherically symmetric space-times. 
Let us consider the following line-element: 



de 



+ p^{T,Odn^ (Bl) 



where A, B, p are continuous, differentiable functions of 

(r, ^) and dfl"^ = dO"^ + sin^ OdcjP' is the metric on the unit 
2— sphere. The action for the scalar fleld propagating in 
the above background is given by 



-9 5^" dfi"^ d^f 



Im 



(B2) 



where we have decomposed <^ in terms of the real spher- 
ical harmonics {Zim{6 , (j))): 



(B3) 



Im 



Following the standard rules, the canonical momenta and 
Hamiltonian of the fleld are given by 



H,. 




(B4) 

AB p\d^^,y {^h) 

H = Him ■ 



The canonical variables (</?;„, H,^) satisfy the Poisson 

brackets 

w,AT,o,u,jT,a} = s{^-a (B6) 

(r, 0, (r, a} - = {n,„^ (r, 0, H,„^ (r, C')} • 

Having obtained the general Hamiltonian, our next step 
is to show that this reduces to the flat space-time Hamil- 
tonian of the scalar field in a fixed Lemaitre time. In the 
time-dependent Lemaitre coordinates [22, 27] the line- 
element is given by (Bl) with 



A{t,0 = 1; b{t,^) = 



; p{t,0 



(B7) 



where r ~ r(T, ^). 

The line-element in Lemaitre coordinates is related to 
that in the time-independent Schwarzschild coordinates, 
viz., 

ds^ = -f(r)dt' + ^ + r'^dV? ■ f{r = rh) = (B8) 



13 



by the following transformation relations [27] : 



t± dr 



fir) 



^ = t+ Idr 

dr 



[l-/(r)]-V2 



(B9) 



The advantage of the Lemaitre coordinate over the 
Schwarzschild coordinate is that (i) the former is not sin- 
gular at the horizon as opposed to the latter, and (ii) 
^(or, r) arc spacc(or, time)-like everywhere while r(or, t) 
is space(or, timc)-like only for r > Vh- 

Substituting the relations (B7) in the general Hamil- 
tonian (B5), we get, 



1 r 



1 



+ 



(BIO) 



where the conjugate variables satisfy the Poisson brackets 
(B6). Note that the scalar field and the Hamiltonian 
depend explicitly on the Lemaitre time. 

Next, choosing the the Lemaitre time (r = Tq = 0), 
the relations (B9) lead to: 



dr 



1 



(Bll) 



If we set d9 = dcf) = 0, then for the fixed Lemaitre time 
To it follows that ds^ = d^'^ /B{to,^) = dr'^, i.e., the 
covariant cut-off is \ds\ = dr. 

Substituting the above relation (Bll) in the Hamilto- 
nian (BIO) we get 



2 70 



dr 



n2 



+ r' {d,^^J^ + l{l + l)cpl 



(B12) 



where the variables (<p,^ , H,^ ) satisfy the relation: 

{^,„(r),n„„(r')} = V^~fir)Sir-r'). (B13) 



Performing the following canonical transformations 



one obtains [28] 



H 



1 r'^ 

pi 



d_ f (fim{r) 
dr 



(B14) 



(B15) 



This is nothing but the Hamiltonian of a free scalar field 

propagating in flat space-time. This is true for any fixed 
value of T, provided the scalar field is traced over either 



the region r G (0,ni] or the region r € [r/i,oo). Note 
that the black-hole singularity can be entirely avoided for 
the latter choice. Now for evaluating time-independent 
quantities such as entropy, it suffices to use the above 
Hamiltonian (the same cannot be said for time-dependent 
quantities). 

The approach here differs from that of Ref. [26] where 
the authors divide the exterior region r > into two by 
introducing an hypothetical spherical surface and obtain 
the entanglement entropy of that surface. In contrast, 
we consider the complete r > rg region and obtain the 
entropy for the BH horizon. We discuss the possible ex- 
tensions in Sec. (VI). 



APPENDIX C: ENTANGLEMENT ENTROPY 
FOR GS AND ES 

For the sake of completeness, we outline the essential 
steps in the computation of entanglement entropy for 
ground state and first excited state. In the following, 
we denote all the quantities, viz., wave function, density 
matrix etc. by the symbol/suffix for GS and by 1 for 
the first ES. 



1. Ground state 

In this case the wave function (10) reduces to (on set- 
ting Ui = 0, for all i): 



{xu...,xr,) = n^(°)exp(-i4/^x?) 



TT 



1/4 



exp 



X fix 



Let us decompose 



n = 



A B 
C 



and define 



B^A-^B 



7 



C-(3, 



(CI) 



(C2) 



(C3) 



where A is an n x n symmetric matrix, i? is an n x (N — n) 
matrix, and C, (3, 7 are all {N — n) x {N — n) symmetric 
matrices. The density matrix (12) reduces to [7]: 



rN-n\ 



A\ 



exp 



(C4) 

In the above, the matrices B and ^ are non-zero if and 
only if the HOs are interacting. Due to the Gaussian 
nature of pf){t; t') (C4), one can make a series of unitary 
transformations : 



V-fV'^ = 7D = diag , /3 = i-o'^VpV^-f-^'^ 
WPW'^ = f3D= diag , v = W^j]^^V , 



(C5) 
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such that it reduces to a product of (N — n), two HO 
{N = 2) density matrices, in each of which one osciUator 
(n = 1) is traced over [7]: 



Po(t;0 = 



11^1 



tN- 



exp 



In the above, the matrices B and /3 are non-zero if and 
only if the HOs are interacting. Due to the Gaussian 
nature of po(t; t') (C4), one can make a series of unitary 
transformations : 

VyV^ = 7d = diag , = 7d'^^V(3V^j-'/^ , 

WpW'^ = pD = diag , V = W'^^]l'^V , (C7) 

such that it reduces to a product of (N — n), two HO 
(A'' = 2) density matrices, in each of which one oscillator 
(n = 1) is traced over [7]: 

N-n 



Po{t;t') 



n '^'^p 



(C8) 

where Vi £ v and Pi £ (3. The corresponding entropy is 
a sum of two HO entropies [7]: 



Pi 



(C9) 



For the total Hamiltonian H, the entanglement entropy 
is 



1=0 



S=^{2l + l)Si=0.3(- 



(ClO) 



where (2^ + 1) is the degeneracy factor that follows from 
the spherical symmetry of the Hamiltonian, R = a{n + 
0.5) is the radius of the hypothetical spherical surface 
the horizon - the DoF inside of which are traced over, and 
Si is the entropy for a given I. Although ideally the upper 
limit Imax should be infinity, for numerical estimation of 
the entropy (for a certain precision) a very large value 
of Imax is assigned in practice. The precision goal Pr is 
set by demanding that the maximum value of I = Imax 
should be such that the percentage change in entropy 



S{1 

max 

)-S{l 



5) 



5) 



X 100 < Pr . 



(Cll) 



For instance, if Pr = 0.01, the numerical error in the 
total entropy estimation is less than 0.01%. 

In the cases of GCS and SS, it has been shown in [14] 
that the expression for the total entropy is the same (up 
to irrelevant multiplicative factors) as that for GS. 



2. First excited state 

In this case, we consider one HO is in the excited state 
while the rest - 1 are in their GS [13, 14]. From 



Eq. (10), we have 

N ,^ ^ 1/4 



X exp 



. (C12) 



In terms of the pure GS wave function (CI) this can be 
written as 

V'i(a;i . . . xn) = \/2 {oFKy,"^^ ipo{xi,..., xn) , (C13) 
where 

= {ai, . . . jUn) (C14) 

are the expansion coefficients and the normalization of 
ipi requires a^a= 1. 

Using Eq.(12) density matrix can be evaluated and is 
given by: 

/n 
JJ dxi [x'^Ax] Vo {xi-, t) Vo [xi-, t') 



Po{t;t') 
(C15) 



where A is a. N x N matrix defined by 

A = U-K]/^aa-K^J'u^^l^^ (C16) 

is an n X n symmetric matrix; As is an n x (N — n) 
matrix; Ac is an (A^ — n) x {N — n) symmetric matrix, 
and K = Tr(AAA~^). The (N - n) x (TV - n) matrices 
A/3 and A^ are given by 



1 



A^ = - 2Ac 

- B^A-'A 

1 



AlA-'B 



b'^a-'AaA-'b) 



-{2AlA-^B - B^A-^AaA-^B) . (C17) 

K 

Af} is symmetric, while A^ is not necessarily symmetric 
due to the presence of the first term on the right hand 
side. 

In general, unlike the GS density matrix po (C4), the 
ES density matrix pi (C15) cannot be factorized into 
(N — n) two HO density matrices. However, po is a 
Gaussian that attenuates virtually to zero beyond its few 
sigma limits. Therefore, if 

£l = tmax-^pimax ^ , ^2 = tmax^l^niax 1 (C18) 

where 

3(iV-n) 



= 

''max 



V2Tr(7 - /?) 



(1,1,...) 



(C19) 
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corresponding to Scr limits of the Gaussian inside po, then 

one may approximate 



1 - 



\-t Apt 



exp 



(C20) 



(C21) 



As such, with the following shift of parameters 

P' = P + A0, j' = j + A^ 
the approximated ES density matrix is also a Gaussian 



pi{t;t') « AC exp 



(C22) 



This can be factorized once again into two HO density 
matrices, and the associated entanglement entropy can 



be computed. For the following set of values: N = 
300, n 100 - 200, o = 10 - 50, o being the num- 
ber of last non- vanishing entries in the vector a"^, i.e., 
= {l/^/o) (0, • • • , 0, 1, • • • ,1), the entropy computa- 
tion is done numerically (using MATLAB) in [14]. The 
precision setting in the computation had been 0.1% and 
the criteria (C18) is satisfied for the above choice of the 
parameters N,n and o. The results show that the ES 
entropy scales as a power of the area. The power, how- 
ever, is always less than unity (for any o > 0) and is 
lesser and lesser, the higher the value of o. The AL is 
thus always violated for the chosen range of values of 
n(= 100 — 200). The other interesting observation made 
in [14] is the shifting of peaks in the variation of the par- 
tial waves {21 + l)Si with I for ES as compared to the 
case for GS. 
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